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' We discuss the BCS-BEC crossover in a degenerate Fermi gas of two hy- 
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I ' perfine states interacting close to a Feshbach resonance. This system has 

^ I quasi-molecular Bosons associated with a Feshbach resonance, and this kind 

■ 

I— 1 1 of coupled Boson- Fermion model has been recently applied to Fermi atomic 

I gases within a mean-field approximation. We show that by including fluctu- 
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\^ ' ation contributions to the free energy similar to that considered by Nozieres 
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I and Schmitt-Rink, the character of the superfluid phase transition continu- 
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I ously changes from the BCS-type to the BEC-type, as the threshold of the 
> ' 

' quasi-molecular band is lowered. In the BEC regime, the quasi-molecules be- 

B' 

1 ' come stable long-lived composite Bosons, and an effective pairing interaction 
■ 

Q I mediated by the molecules causes pre-formed Cooper-pairs to also be stable 

O ■ 

^ I even above the transition temperature Tc. The superfluid phase transition is 

^ ' shown to be usefully interpreted as the Bose condensation of these two kinds 

■ - - ' of Bosons. The maximum Tc is given by (in terms of the Fermi temperature 

Tp) Tc = 0.218Tf (uniform gas) and Tc = 0.518Tf (trapped gas in a harmonic 
potential). 
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The achievement of a superfluid phase transition in a trapped gas of Fermi atoms is one 
of the most exciting challenges in the current study of ultracold atoms. At present, several 
groups have produced degenerate Fermi gases down to around 0.2Tf (where Tp is the Fermi 
temperature). |[T|] In the past few years, many theorists have discussed the properties of a 
superfluid phase based on the BCS concept of Cooper-pair formation. More recently, 
this approach has been extended to deal with a Feshbach resonance. A new strongly 
attractive interaction can arise between Fermions mediated by the quasi-molecular Boson 
associated with a Feshbach resonance. However, one must be careful in applying BCS theory 
to a Fermi gas when the pairing interaction is very strong. In this case, fluctuations in the 
two-particle Cooper channel strongly suppress the transition temperature Tbcs predicted by 
weak-coupling BCS theory. Moreover, in this strong coupling regime, one has two "kinds" 
of Bosons even in the normal phase above T^, i.e., molecules associated with the Feshbach 
resonance and pre-formed Cooper-pairs. Thus the superfluid phase transition is expected to 
make a crossover to Bose-Einstein condensation (BEC) from the BCS-like one, which 
is the subject of the present letter. 

Our starting model is identical to that developed in Refs. but we extend their 
mean-field analysis to include fluctuations in the particle-particle channel. For simplicity, 
we discuss a uniform Fermi gas. Our work clarifies how the character of the superfiuid 
phase transition gradually changes from the BCS-like transition to BEC-like transition as 
the threshold energy for the formation of quasi-molecular Bosons is lowered. 

We consider a uniform gas of Fermi atoms with two hyperfine states which interact close 
to a Feshbach resonance. As usual, we describe these states as pseudo-spins a =t, i, with 
N-^ = Ni = N/2. Our model system is described by the Hamiltonian 

per ppl 

+ Y^i^El + 2u)h% + ^7 Ei^IW2TC-p+,/2i + h.c]. (1) 
q p,q 

Here Cp^ and represent the annihilation operators of a Fermi atom with the kinetic energy 
Ep = /2m and a quasi-molecular Boson with the energy spectrum E^ + 2v = /2M + 2z/, 
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respectively. The threshold energy of the composite Bose particle energy band is denoted by 
2v and is the crucial adjustable parameter in this model. The attractive pairing interaction 
—U < originates from non-resonant processes, while (7 is a coupling constant associated 
with the Feshbach resonance. Since the 6-Boson is made of two Fermions, we have M = 2m 
and impose the conservation of the total number of bare Fermi atoms: {N) = {J^pa Cpo-Cpo-) + 
2{J2qblbq) . This crucial constraint is incorporated in the Hamiltonian when we replace H 
by H — fiN. The resulting "grand-canonical" Hamiltonian has the same form as (|l]), except 
that now Sp Ep — /ip and — jj,^, with /ip = /Ub/2 = /i. A coupled Fermion-Boson 

Hamiltonian similar to (0) has been extensively discussed in connection with the high-Tc 
cuprates, usually setting f/ = 0. 

The superfluid phase transition temperature can be conveniently determined by the 
Thouless criterion. This describes the instability of the normal phase of Fermions with an 
attractive interaction associated with the formation of Cooper-pairs of atoms with opposite 
spins. We first calculate the four-point vertex function r(j9+,p_;p'_^,p'_) within a gen- 
eralized t-matrix approximation in terms of U and g (as shown in Fig. 1). One sees that 
this involves an infinite sum of ladder diagrams involving U and particle-particle bubble 
diagrams involving g. The equation for the particle-particle vertex F is given by 

r(p+,p_;pV,P-) = (/^o(g)-f/) 

+ {g^D,{q) - U)T Goipl)Goipl)T{pl,p'L;p'^,p'_), (2) 

P 

where p+ = (p+q/2, ium+ii^n), P- = (p— q/2, iuJm) and q = (q, iun), with ium and iz/„ being 
the Fermion and the Boson Matsubara frequencies, respectively; Gq ^(p, iciJm) = icom — £p + fJ' 
and DQ^{q,ih'n) = i^n — -E^ — 2z/ + 2/i are, respectively, the bare one-particle Fermion and 
the Boson thermal Green functions. The solution of (0) is 

r q,«z/„ = — — — . — — . (3) 

1-[U - g^Do{(i, «z/„)JH(q, zz/„) 

Here H(q, iz/„) is the correlation function of the Cooper-pair field operator of two atoms with 

total momentum q defined by A(q) = Z]p Cp^q/2T'^-p+9/2i' absence of U and g: 
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n(q, 11.^) = y: (4) 

p ^p+q/2 + ^p-q/2 — 2/i — iVn 

where f{e) is the Fermi distribution function. n(q, iz/„) is often referred to as the particle- 
particle propagator describing superfiuid (Cooper-pair) fluctuations. 

The superfiuid phase transition occurs when the particle-particle vertex function in 
develops a pole at q = z/„ = 0. This is easily seen to correspond to the following equation 
for Tc. 

In (H), (7^/(2z/ — 2jj) (= —(7^/^0(0)) can be interpreted as the additional pairing interaction 
mediated by a molecular Boson. This interaction becomes very large as 2/i 2v. Eq. @ 
defines the temperature Tc at which an instability arises in the normal phase of a Fermi 
gas described by (|lD due to formation of bound states with zero center of mass momentum 
(q = 0) and energy 2/i. 

In Refs. [§,0, the model in (|l]) is solved in the superfiuid phase within a simple mean- 
field approximation in which the 6-Bosons are replaced by a Bose condensate in the q = 
state, described by (&,=o) = 0m- Within the BCS-pairing approximation to the two-particle 
interaction, the Hamiltonian H — fiN is then easily diagonalized. One finds the standard 
self-consistent BCS equation for the effective gap function defined by A = —Up—g(j)m, where 
p = (Sfe CfctC-fcj^) describes the Cooper-pair Bosons, with the effective attractive interaction 
Ucs = U + /{2v — 2/i). In this mean-field BCS theory below Tc, the quasiparticle energies 
are given by Ef, = \J {e^ — /i)^ + | Ap. Within this mean-field theory, the condition that the 
gap function A vanishes can be shown to lead precisely to (§) as given by the Thouless 
criterion. Thus (|) gives a very clear picture of the physics discussed in Refs. [|,|]. 

The chemical potential /i for our model is determined from the condition that the total 
number of particles (Fermion+2x Boson) in the interacting system is equal to the number of 
bare Fermions A^. As well known, this simply gives /i = in the weak-coupling BCS theory 
if we neglect the temperature dependence of /i (ep is the Fermi energy.). In the mean-field 
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approximation of Refs. @J^, one has simply = A^p + 2|0mP below Tc. At Tc, we have 
(pm = 0. However, /i(Tc) in (|^) can be different from ep simply because the temperature 
dependence of /i(T) in a non-interacting Fermi gas is significant when T is comparable 
to the Fermi temperature Tp. It is clear that in Refs. the transition temperature 

corresponds to breaking up bound states with q = 0, rather than exciting these molecules 
into states of finite q. However, /i can deviate from ep in a more fundamental way when the 
quasi- molecules with q 7^ 0, superfluid fluctuations, and the pre- formed Cooper-pairs are all 
properly included. [l6|,[7|,p!0[| This is the goal of the present letter. 

The total number of Fermi atoms N{fi, T) can be obtained from the thermodynamic 
potential Q using the identity A^ = —dQ/d^. The effect of the fluctuations in the two- 
particle vertex function on the potential Q is shown in Fig. 2. [^] In Fig. 2, (a) is the 
same as the fluctuation free energy diagrams which were studied in Ref. 0. In contrast, the 
diagrams in (b) originate from the Feshbach coupling of the 6-Bosons and Fermions in (P, 
and are new in the present analysis. One finds [q = (q, iz/^)] 

(6) 

q ^1^ 

Here A^^ = 2Y.pf{ej,- /i) and = E5^b(^° + 2v - 2/i), where n^{E) is the Bose dis- 
tribution function. For g = U = 0, the 6-Bosons in (1) start to Bose condensate when 
2z/ < 2ep (see Fig. 3(a)). The last term in (|^) comes from adding the two kinds of thermo- 
dynamic fluctuations shown in Fig. 2. Note that the fluctuations in (|^) are the same as the 
denominator of the vertex function in (^. We leave out the convergence factor e^^'^" in the 
subsequent frequency sums. 

We must solve (|^) together with (^ self-consistently. The results we obtain for /i(Tc) 
and Tc as a function of the threshold energy 2z/ are shown in Figs. 3(a) and 3(b). However, 
before discussing these numerical results, it is useful to first give a physical interpretation 
of the contributions in (^ in terms of the effect of U and g on two kinds of Bose molecules 
(two-particle bound states). We proceed by using the identity = —TJ2q,un Doi^qjii'n) to 
rewrite @ as 

5 



A^ 



A^^ 



r) 1 1 



[U-g^D,{q)]Ii{q) 



N = N^- 2TY^b{q) - -— ln[l - U^M^io) 

q q C/^F L 



(7) 



Here D~\q) = zz/„ - E° - 2, 



^n(g) is a renormalized 6-Boson Green function, with 



n(g) = n(g)/[l — f/n(g)]. The second term 2A^b in (0) can be interpreted as the contribution 
of the 6-Bosons as affected by the self-energy = -g^U{q). The third term 2Nc in 
has the same form as the contribution of fluctuations studied by Nozieres and Schmitt-Rink, 
except that now the effective attractive interaction Ues{q) = U — g^Do{q) depends on 
energy as well as momentum. Thus this third term 2Nc may be viewed as the effect of 
Cooper-pair fluctuations on A^. We note that the mean-field approximation discussed in 
Refs. is reproduced by (|^) by ignoring 2Nc, as well as setting g^Il{q) = in and 
only considering the q = &-Boson (i.e., setting = 0). In the strongly coupled Fermion- 
Boson problem we are dealing with, it should be emphasized that our interpretation of the 
separate terms in as 2A^b and 2Nc is only meant as a way of understanding the essential 
physics. 

In understanding the character of the renormalized 6-Bosons described by A^b in 
it is useful to transform the Matsubara frequency summation into a frequency integration, 
Nb = —^J2qI^oo^^''^B{z)I'mD{q,ih'n ^ z + i6). Here we take the principal value in the 
^-integration. One sees that if the 6-molecule decays into two Fermi atoms in the presence of 
the Feshbach resonance, it has a finite lifetime as given by the inverse of the imaginary part of 
the self-energy —g^H in Z)(q, z+i6). However, when the chemical potential becomes negative 
by lowering the threshold energy 2i/, a new situation arises: Since it can be shown from 
that /mH(q, z + i6) is proportional to the step function Q{z + 2/i — E^), the renormalized 
6-molecules do not decay if their energies are smaller than E^ — 2fi. The energy {uq) of a 
stable molecule corresponding to the pole of D is given by (for E^ — 2fi > 0) 



Such stable long-lived molecules appear when the "effective" threshold 2z/ = 2z/ — (7^H(q, z) 



u;, = (E;-2/i) + [2z/-/n(q,cu,)]. 



(8) 
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becomes negative, since decay into two Fermi atoms is forbidden. Around 2; = q = 0, 
we may use the static approximation for this renormahzed threshold, 2z> ~ 2z/ — g^ll{0, 0). 
Using (^, one can show that 11(0,0) > 0, and hence 11 is found to promote the generation 
of stable molecules by decreasing the value of 29. 

N-Q in (|^) is thus found to consist of two kinds of Bosons, A^b = + N^^^ , i.e., 

stable molecules N^~^ with infinite lifetime and quasi-molecules A^g^° which can decay 
into two Fermi atoms (Strictly speaking, N^^^ also includes a contribution from scattering 
states.). Extracting the contribution of the poles describing these stable molecules, we obtain 
A^B~° = X]q°^'^^ Z{q)nB{uJq), where Z{q)^'^ = 1 + (^^^n^o^ ^ggf^jj^ji^gg j^iass renormalization. 

A similar discussion can be used to understand the third term 2Nc defined in (|^), de- 
scribing the Cooper-pair fluctuations. When one writes the non-resonant s-wave interaction 
U in terms of the s-wave scattering length U = AiraN/m, one finds U/e-p = 0.85(pFa), 
where pp is the Fermi momentum. Since a dilute Fermi gas is characterized by ppa <^ 1, 
we must limit ourselves to U/ep <^ 1. In this case, the non-resonant attractive interaction 
U in (0) cannot give rise to the pre- formed Cooper-pairs (in a three-dimensional system). 
P|7|JT0[] However, when 2/i approaches 2z/, the interaction Ues{q) mediated by the molecules 
becomes very strong, so that stable pre-formed Cooper-pairs may appear. The energy of 
these poles is found to be precisely the same as that given in (||). Thus we can also divide 
A^c iiito contributions from (stable) pre-formed Cooper-pairs (= Nq~^) and from scattering 
states (= A"^^), with N^=^ = ^ ^2gf^Z(q)nB(cu,). We see that if ^ = or if we 

ignore fl, there are no stable pre-formed Cooper-pairs and Z{q) = 1 in Nq~^. 

These stable Bosons give a very clear physical picture of BCS-BEC crossover present in 
the model Hamitonian in (|l|). When the threshold 2z/ is much larger than ^f, the Fermi 
states are dominant. In this regime, a BCS-like phase transition is found. Indeed, since 
yU ~ £:f > in this regime, stable Bosons are absent because Jmn(q, z) ^ for z > 0. The 
only exception is for z = (since /mn(q, 0) = 0). However, the condition (2z> = 2fi) that 
a stable Boson with q = appears at cj, = given by (^) is easily shown to reduce to the 
equation for Tc in (|^). Namely, as well known in weak-coupling BCS theory, the formation 
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of stable Cooper-pair Bosons and the phase transition occur at the same temperature. 

In this hmit, no stable long-lived Bosons with q 7^ exist above this transition temperature. 

In the opposite limit z/ ^ 0, the Fermi states are almost empty. In this regime, the 
phase transition is expected to be BEC-like. Indeed, since the chemical potential should 
be negative for z/ <^ 0, stable Bosons can appear even above T^, as discussed above. The 
phase transition of the gas of these stable Bosons occurs as usual when the energy of the 
Boson with q = measured from the chemical potential 2/i reaches zero. Eq. @) with this 
condition {ojq = q = 0) again reproduces (H), i.e., 2/i = 2i> = 2z/ — (7^11(0,0) is equivalent 
to 1 = f/cfr(0)n(0, 0). In the limit ^ 0, the problem has effectively been reduced to the 
BEC transition in a non-interacting gas of A^/2-Bosons with mass M = 2m, with no free 
Fermions. In this limit, since the superfluid Cooper-pair fluctuations described by 11 (g) are 
absent, @ [or (||) with cj, = q = 0] simply gives 2/i = 2i>, the required condition for a 
Bose condensate to appear in N^~^ and Nq~'^. Thus, although was originally obtained 
as the condition for the superfluid phase transition in a Fermi gas due to the formation of 
two-particle bound states, it also describes BEC in a gas mixture of stable 6-molecules and 
pre-formed Cooper-pairs {N^~^ + Nq~^). 

We now discuss our numerical solution of the coupled equations and (|^), which can 
be understood in terms of preceding analysis. As usual, a cutoff in calculating n(q, iz/„) is 
needed to remove an (unphysical) ultraviolet divergence. We simply introduce the factor 
exp[— (er/cUc)^] in the energy integration, with = 26^. We give energies in units of the bare 
Fermi energy ep and particle number in units of the total number of bare Fermi atoms. 
Fig. 3 (a) shows the dependence of Tc on the threshold parameter 2z/. As u is lowered, Tc 
deviates from the BCS limit. At the same time, at Tc the character of the particles involved 
gradually changes from Fermions to composite Bosons. This is shown in Fig. 3(c), as 2/i 
approaches 2z/ in the BEC limit (Fig. 3(b)). We also see in Fig 3(c) that the character 
of the Bosons also changes from a "Feshbach-resonance" Boson (quasi-molecule A^b^*^) to a 
mixture of stable 6-molecules and pre-formed Cooper-pairs (A^b~° + Nq~^) as 2z/ decreases. 
I 111 In the limit z/ <^ 0, we see that the pre-formed Cooper-pairs disappear, and system 
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reduces to a gas of 6-Bosons. 

In the BEC limit, since the Fermions, pre- formed Cooper-pairs and superfluid fluc- 
tuations are all absent, Tc is simply determined by the equation = 2N^. Using 
N = (2mTF)=^/V3v^^ we obtain Tc = 2Tf/[QV^C{^/'^)?^^ = 0.218Tf. This value of 
Tc gives the maximum phase transition temperature within our model [U/ep < 0.3). 

The BCS-BEC crossover induced by changing the strength of the pairing interaction has 
been investigated as a model for high-Tc superconductors. 0Jl3l In the present case, when 
we plot Tc with respect to the effective pairing interaction, Ues = U + g'^/{2v — 2/i), we 
obtain the inset in Fig. 3(b). This is analogous to results in Ref. [0]. Since 2/i approaches 
2v as V decreases (see Fig. 3(b)), the strong coupling regime in the inset corresponds to the 
case of small v in Fig. 3(a). 

For trapped Fermions, one expects that the maximum Tc will again be given by the BEC 
transition for A^/2-Bosons but now in a parabolic trap. This gives Tc = Tf(1/6C(3))^/^ = 
0.518Tf, where Tf = {?,Nfl'^ujQ. 

To summarize, we have generalized the analysis in Refs. to show that the superfluid 
phase transition in a uniform gas of Fermi atoms with a Feshbach resonance can change from 
the BCS-like to BEC-like when we include the fluctuations associated with the Cooper-pair 
(or two-particle) channel. We have also introduced a physical interpretation of our results 
in terms of two kinds of Bose molecules. In contrast with superconductors, the BCS-BEC 
crossover in atomic Fermi gases appears to be accessible by lowering the threshold energy 
2v of the quasi-molecular Bosons associated with the Feshbach resonance. Generalization of 
these results to below Tc and for a trapped Fermi gas within the local density approximation 
(LDA) will be discussed elsewhere. 
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FIGURES 

FIG. 1. Generalized t-matrix approximation for the two-particle vertex function F. The solid 

(dashed) line represents Gq {Dq), the dotted line describes —U, and the vertex g is the coupling 
parameter defined in (1). 

FIG. 2. Thermodynamic potential 0, fluctuations associated with the two-particle Cooper-pair 
channel. Infinite series of these kinds of diagrams are summed. 

FIG. 3. Self-consistent solution for (a) Tc and (b) fi{Tc) as a function ofv,U = O.Sep, 9 = 0.6£f, 
and a cutoff uJc = 2eF- In panel (a), for comparison, we show the BEG case = ?7 = 0), and BCS 
solution [4] of (5) with /u(T) being the chemical potential of a non-interacting Fermi gas. Panel (c) 
shows the number of Fermions iVp, stable 6-molecules N^~^, quasi-molecules (finite lifetime) N^^^, 
pre-formed Cooper-pairs N^~^ , and the scattering part Nq, all at Tc, as a function of the Feshbach 
threshold parameter 2z/. In the case shown, stable long-lived Bosons only exist for v <, 0.22£f. 
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Fig. 1 Ohashi and Griffin 




Fig. 2 Ohashi and Griffin 




Fig. 3 Ohashi and Griffin 



